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Abstract: Gravity in five-dimensional braneworld backgrounds often exhibits problematic 
features, including kinetic ghosts, strong coupling, and the vDVZ discontinuity. These 
problems are an obstacle to producing and analyzing braneworld models with interesting 
and potentially observable modifications of 4d gravity. We examine these problems in a 
general AdS^/ AdS^ setup with two branes and localized curvature from arbitrary brane 
kinetic terms. We use the interval approach and an explicit "straight" gauge- fixing. We 
compute the complete quadratic gauge-fixed effective 4d action, as well as the leading 
cubic order corrections. We compute the exact Green's function for gravity as seen on 
the brane. In the full parameter space, we exhibit the regions which avoid kinetic ghosts 
and tachyons. We give a general formula for the strong coupling scale, i.e. the energy 
scale at which the linearized treatment of gravity breaks down, for relevant regions of the 
parameter space. We show how the vDVZ discontinuity can be naturally but nontrivially 
avoided by ultralight graviton modes. We present a direct comparison of warping versus 
localized curvature in terms of their effects on graviton mode couplings. We exhibit the 
first example of DGP-like crossover behavior in a general warped setup. 
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1. Introduction 



Gravity in five-dimensional braneworld backgrounds exhibits features whicli, from a 4d 
point of view, are both novel and surprising. Examples include: 

• Kaluza-Klein (KK) graviton modes whose couplings to brane matter are only TeV 
suppressed Q, and thus potentially detectable as resonances at the LHC 0-0] • 

• Extra 4d scalar modes (radions) whose couplings to brane matter are TeV suppressed 



but otherwise like those of gravity |]5|-[12]. 



Continuum KK graviton modes which are potentially detectable at the LHC [13|,|14|. 



KK graviton modes whose couplings to brane matter are suppressed by warping 

KK graviton modes whose couplings to brane matter are suppressed by localized 
curvature {i.e. brane kinetic terms for gravity) [^]. 



4d effective theories with no massless graviton |16| ,[17| 



4d effective theories with a massless graviton and an additional ultralight KK gravi- 
ton, whose mass goes to zero in a well-defined parametric limit |jl^ 



A crossover scale, such that KK gravitons with mass below this scale have unsup- 
pessed couplings to brane matter, but heavier KK modes have suppressed couplings. 
This leads to modifications of 4d gravity which appear in the infrared rather than the 
ultraviolet [17|,[20|-|22]. Such "DGP" scenarios could have relevance to cosmology 



What is truely remarkable is that these novel features of gravity do not require the 
assumption of any new physics beyond a single extra dimension and the existence of co- 
dimension one branes. The above results are derived in low energy linearized effective 
descriptions which require no assumptions about quantum gravity, string theory, or the 
existence of exotic matter. 

Before we celebrate too loudly, however, we should confront a number of troubling 
issues common to most or all braneworld gravity models: 

• Are these low energy effective descriptions really just special cases of 5d General 
Relativity? This is not obvious since most models invoke orbifold backgrounds, and 
since they also ignore the problem of ultraviolet matching to e.g. the microscopic 
features of the branes. 

• What are the physical degrees of freedom in these various models? This is a difficult 
question to answer rigorously, since the usual gauge- fixing choices {e.g. harmonic 
gauge) are not suitable for setups with more than one brane. 
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Are these models stable? Models with massless radions are at best marginally stable. 
Some simple setups [26|-||2^ have radions which are kinetic ghosts, i.e. their kinetic 



terms have the wrong sign. Kinetic ghosts indicate an instability in the model. As 
we will see in this paper, graviton modes can also be kinetic ghosts, and some kinetic 
ghosts are also tachyons. Although kinetic ghosts may be useful for some purposes 
their presence cannot safely be ignored. 



Under what conditions do ultralight graviton modes mimic 4d gravity, avoiding the 
famous vDVZ discontinuity |34]? 



At what scale does the low energy linearized approximation break down? Perturba- 
tion theory involves more input parameters than just the 5d gravitational coupling 
M; it would not be surprising if in certain limits of these additional parameters the 
low energy effective theory breaks down at a scale much lower than M. Indeed this is 
precisely what happens in the simplest models containing ultralight graviton modes 
or crossover behavior, where strong coupling sets in at a scale parametrically much 
smaller than M. 



In a recent investigation [35| three of us resolved the first two issues listed above. 
We recast the standard braneworld models in the "interval picture", where orbifolding is 
replaced by intervals with boundaries. We showed that braneworld gravity has a well- 
defined action principle only if we extend General Relativity to include "brane-boundary 
equations" which supplement the usual bulk Einstein equations. We also showed how to 
rigorously extract the physical degrees of freedom, by introducing a class of "straight" 
gauges suitable for braneworld analysis. 

The purpose of this paper is to address the remaining three issues listed above, and 
thus to revamp the promising field of braneworld gravity. In section 2 we derive the general 
4d quadratic effective action for a general setup with two branes, including warping and 
localized curvature. This action is gauge-fixed to just the physical degrees of freedom, using 
a straight gauge. In section 3 we use this effective action to demonstrate the presence 
or absence of kinetic ghosts and/or tachyons according to various choices of the input 
parameters. This analysis is a generalization and improvement of earlier attempts p6|. 



1 37]. We find that both the radion and the graviton zero mode can be kinetic ghosts in 
certain regions of the parameter space. 

In section 4 we address the question of strong coupling. By computing the relevant part 
of the 4d cubic effective action, we give a general formula to determine the strong coupling 
scale for the radion. We show that the strong coupling scale becomes small in a DGP-like 
limit. In section 5 we give an exact expression for the straight gauge graviton Green's 
function on the brane. This allows us to show how the vDVZ discontinuity is avoided in 
warped models with ultralight graviton modes. In a DGP-like limit, our straight gauge 
graviton Green's function does not have any diverging tensor structures such as those that 
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appear in DGP; the potential breakdown of linearized gravity in this limit is entirely due 
to the radion. 

Finally, in section 6 we use the exact Green's function to study the couplings of KK 
gravitons to brane matter. In the special case of the Karch-Randall setup, we regain 
the recent results of Kaloper and Sorbo ||3^. We focus on models with an infinite extra 
dimension and no massless graviton. Localized curvature and warping both have the effect 
of making an ultralight graviton from the first massive KK mode. We compare these two 
effects in our models, and show that the warping effect is more efficient than localized 
curvature in creating a simulacrum of 4d gravity. 

In section 6 we exhibit explict crossover behavior for models in our general setup. In 
these models the couplings of the KK graviton modes to brane matter are unsuppressed 
up to a mass scale 1/rc, and are highly suppressed for modes heavier than this scale. This 
behavior is the most interesting phenomenological feature of DGP braneworld gravity; here 
we can study it in a general warped framework. 



2. The 4d quadratic effective action 



In [35 1, a brane world gravity theory with the action 

S = J d^x (^J^ dy + j dy^ ^^{2M^R- 

+ Z] / d^x^f^){2Mfn^'^ -Vi} + AM^ (f K. (2.1) 

• J y=yi J dM 

was analyzed. This action represents a general warped gravity setup with codimension one 
branes, written in the interval picture. In (|2.1| ) M is the 5d Planck scale, A = —2AM'^k'^ 
is the bulk cosmological constant giving a bulk curvature k, the Mj are the coefficients of 
brane-localized curvatures iZ^'^\ the Vi are brane tensions, and K is the extrinsic curvature 
of the Gibbons-Hawking boundary term. 

Upon linearization, Gmn = G^,in + ^mn, the AdS^/AdS^ background solution is 

= ('l" °) , (2.2) 



where 



with 



_ <y? „ .2 3) 



, , cosh k(v — yo) ^ , 

V^Liy = diag(— 1, 1, 1, 1) and = rj^i^x^x'^. Hereafter, we will omit the superscript ^. 
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Figure 1: The unshaded area is where — 1 < T,- < 1. On the boundary, the curved hues are 



described by Wi 



and the straight hues are Wi — ±12. 



It is extremely convenient to trade Mi and Vi for dimensionless input parameters Vi 
and Wi defined by 

Vi = kMf/M^ , Wi = Vi/2kM^ . (2.5) 

Tlie brane separation L and the warping parameter yQ are then determined from these 
inputs by the relations 

T^ = ^ + "^{l-T^), (2.6) 

or 



with 



Tt = -{-l±\ll + -wm + vf), (2.7) 



To = tanh/cyo > = tanhA;(L — uq) . (2-8) 

The AdS4 inverse radius of curvature H is given by 

Solving— 1 < < 1, Figure H shows the range of input parameters which gives 74(iS'5/A(iS'4 
backgrounds. 

Using this background solution, we eliminate gauge degrees of freedom and find the 
following 4d physical degrees of freedom: 
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1. At the massive level, we have a KK tower of massive spin-2 particles with 5 degrees 
of freedom each. We can decompose the graviton explicitly into this tower of massive 
modes: 

V = E^S =Y.y^''\y)B^K^) , (2.10) 

q q 

where q labels the mass, q = m?/H^, and 

= -P(_i+^9+4^)/2(^) - -7^Q{-l+^9+^)/2(^) ; (2-11) 



here z = tanh k{y — yo) and the P's and Q's are associated Legendre functions. The 
mass spectrum of modes is determined by solving the determinant equation, 

aobi - aibo = , (2.12) 

with 



ao 



{voq{l - Ti) + (3 + ^9T%)ro}p-2^^g^)/2(-^o) 
+(-5 + V'oTi^) ^(l'3+V9Ti5)/2(-^o) 



'1 - {v,q^l - Pl_\^^y,{-To) - 2P-\^^^/^(-ro)} , 
6o = g^r^{^o^r^Q2_^^^^^^(_j.^) + 2Q;_^^^)/,(-To)},(2.13) 



2. At the massless level, there is a massless spin-2 particle (graviton) /3^iy with 2 degrees 
of freedom and a 4d scalar (radion) '4^{x). Thus h^n is 

V = Pi^^ + a^:i^i(y)V^V^V + 9n2^y2{y)i^ . (2.14) 
The y-dependence of is determined to be 

= a\y)B^,{x) . (2.15) 
The remaining y-dependence is in the 3^j's, which are given by 

yi{y) = c{l-zf + dz-^, y2{y) = -H''c{l-zf -H''dz + -T, (2.16) 
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where 



k J^jL) - ■F(O) 
k aoT{L) — Qi^(O) 



with 



(2.17) 
(2.18) 

"i = '■ TTT^^ , (2-19) 

l + kXiTi ^ ^ 

where we have borrowed a notation from |3^: = — = +1- The above formulas 
show a residual gauge freedom parametrized by a single real function F{y), such that 
r{y) = F{y) and d'{y)=:F/a\y). 



2(1 - OiTi, 



ao - ai 



ikXi{l 



Using the background solutions, we can expand the action ( |2.lD up to the second order 
in an arbitrary metric fluctuation h^'^^ . The bulk part becomes 

cL 



the brane part is, with 4d total divergence terms dropped. 



(2.20) 



-g 



h 



2a? 8 



(2.21) 



y=yi 



and the extrinsic curvature part turns into 



(fx 



— /8a' ~, 4a' ~ 4a', 
-g\ h /i H h 

\ a a a 



a M „ o 



y=L 
y=0 ' 



(2.22) 



2a' r,. Ir,^ 3a'. 

— /ia44 n n,44 H 

a 2 a 

A tilde indicates that the corresponding entity is a 4d quantity constructed with g^^. Note 
that there are 5d-total derivative terms in the bulk part. Due to the finiteness of the 5th 
dimension, they do not vanish identically but make contribution to the brane-boundary 
part of the action. 

Expanding further the bulk part, with the help of 

„/ 



-T, 



(2.23) 
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and imposing the partial gauge choice /i^4 = 0, (|2.20| )+( |2.2lD +( |2.22D becomes 



S 



2AP 



8k 



1 _ 9^/2 



h 



12 



'^""^ 2a2 ' 4 ■ 4a 

H /i44 H 5 '^44 — —hh44 —h n44 H n/l44 



2/ 



2a 



^ ~ 2kTih\4 + kTihh44 



y=yi 



(2.24) 



Note that all the terms linear in h get cancelled, as they should be. 

To get a simpler form, we remove y-derivatives on the fields whenever it's possible. 
For example, 



-q — h 
^4a 



2/ 



4(1 



4(1 



. 3a'^ + /c^a^ ~2 



(Sa^a'' + a\") 



4a2 



(2.25) 



where the first term contributes to the brane-boundary part of the action. This way, we 
eliminate /i^', and hh'44-tevms from the bulk part of the action, to get 



2M3 



2 (Tx dyV^[ - 8r + 



- 2a' 
2^ 



12 



-\h^""h^J + ^ - ^W2/i44 + \h44y^Nuh^''' 



3a-_^2 



3a' r/, 3i/2 
"'^■^^ ~ '2a ~^ ~2a? 
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^ 3 ^ 2a2 



V + 



4a2 



Plugging ( 2.1C| ) and (|2.14| ), we finally obtain 



S 



S 



2M3 2M3 



A 2M3 



2M3 



massless 



(2.26) 



(2.27) 
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where 

S 



2M3 



XV -9 



4XiH^ ^ Ui 



y=yi 



(2.28) 



2M3 



<J>0 i 



d'x^g[-B^^^^'v'Bj^J + —B^i)i^'B^} - -^B^^^f^'B^^j) , (2.29) 



and 



2M3 



massless 



i 

+2 1 d^x£dyV^g{ - ^yi - IhW,' - \Fy,)i,V,i, 



= C(0) j d'^x^g (^5^^^ V^fi^. + ^B'^'B^,) + J d'x^g i^A^ , (2.30) 



with 



A; 



k 



= ^ (2 / J 3^^^^^^- + E[-^(i - -')3^^'^'].=.0 ' 



m 



(?)2 ^ ^2. 



J —Tn 1 — 



J%y(^^'dz + z^W-z')y^'^'>\=y/ 



y=yi 



(2.31) 
(2.32) 

(2.33) 
(2.34) 



A hatted entity is defined using the metric without a warp factor. If we include a source 
term in the action, its couphng to a specific graviton mode wiU show up hke 

= 2M34«) J d\^ ^B^i^f'^V^Bj^J + ■■■ + J d'^xdy^ h^^^T'"'' 

[ d'^x^g -mf^'ij''B^$ + ■■■+[ d'^xdy^g ^f^L^t^T'^- , (2.35) 
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and then we can read off the gravitational couphng constant. In particular, on the 0-brane, 

1 y5)(o) 



(2.36) 



3. Avoiding ghosts and tachyons 
3.1 ghostbusting 

It has been known for some time that brane setups of the type that we are considering are 
sometimes afflicted with unphysical features. Kinetic ghosts, by which we mean wrong-sign 
kinetic terms for physical modes in the 4d effective action, are indicative of an instability, 
similar to the case of their cousin the tachyon. Regions of our input parameter space 
which produce kinetic ghosts are certainly to be avoided if we are interested in static brane 
configurations. A kinetic ghost radion can occur for setups with e.g., vq = vl = and 
wq + wl < 0. Intuitively we also expect a kinetic ghost graviton to occur in cases where 
vq and vl become too negative, i.e. we have too much wrong-sign localized curvature. As 
we will see, the full story is quite complicated. 

The boundary between a region of the input parameter space which has a kinetic ghost, 
and a region which does not, defines a class of models where the coefficient of the kinetic 
term of a physical mode is vanishing. After a canonical rescaling of the field, this implies 
strong coupling once we go beyond the linearized theory. Such regions of strong coupling 
are to be avoided if we want the 4d effective low energy theory to be valid up to energy 
scales approaching k. 

In this section we will map out the input parameter space and identify the region which 
avoids both kinetic ghosts and strong coupling. 

We have already noted that the massless mode of the graviton may be a ghost. When 
we choose {vq, wq) such that Tq > 0, Figure § and ^ show how Cg^^ = kcf'^ /{I — Tq ) varies 
as a function of {vl,wl). Cg^^ is zero along each line shown in Figure |2| and ^ positive 
above it and negative below it. As Tq approaches 1 either by vq +oo or by wq — > +12, 
Cg^^ = line moves to the left. 

We can play a similar game with the coefficient, C^, of the radion. Equation ( |2.34| ) 
holds irrespective of gauge choice of F{y). Then, in a generic case where oq — ai ^ 0, for 
simplicity we can choose 

•^(y) = x(^(aL -ao) +ao) , (3.1) 

where x is a remaining real gauge parameter that we leave arbitrary as a check of general 
covariance for physical results. With this choice we get 

c = — EJ2' d = 0. (3.2) 
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Figure 2: The Cg^'' = on (w^, WL)-plane when vq — 0.3 and wq — 10, 5, —5, —10 for the hnes from 
left to right respectively. 




Figure 3: The Cg°^ = on (wl, WL)-plane when wq = 5 and vq — 1.5,0.5,-0.5,-1 for the lines 
from left to right respectively. 



Then, (2.34) becomes 



k 



, - , + Viz)a, ■ „ 



z 1 — z 
ai + 2- ) + 



1 + z 



1-z 
— ; — Vi 



l + z 



2k 



1 + ViTi 



(3.3) 
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Figure 4: The = on (w^, WL)-plane when vq — 0.3 and wq — 10, 5, —5, —10 for the hnes from 
left to right respectively. 




Figure 5: The = on (w^, WL)-plane when wq — 5 and vq = 1.5, 0.5, —0.5, —1 for the lines from 
left to right respectively. 

Figures § and || show that, if we choose {vq, wq) such that To > 0, then = 2kC^/'ix^H'^ 
is positive above each line. Also, the closer Tq is to +1, the more convex to the left = 
line gets. 
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3.2 tachyons 



Equation (2.12) determines the KK mass spectrum of the graviton. For some range of 
values of the input parameters, it can have a zero at negative g, i.e. negative mass- 
squared. For example, if we choose vq = 5,wo = 3,vl = —7,wl = 9, then ( 2.12| ) has a 



solution q = —0.94. 

We will refer to such solutions as tachyons. The tachyons with —9/4 < q < obey the 
Breitenlohner-Freedman bound |^^, and as expected we will find (see section 5.1) that the 
Euclidean Green's function for such modes has the same exponentially damped asymptotics 
as for ordinary massive modes. The tachyons with q < —9/4 have oscillatory asymptotics, 
similar to tachyons in flat space. It turns out that some of the tachyonic solutions also 
have wrong-sign kinetic terms. 

We can expose the tachyonic solutions by analyzing the quantity 

q{l-Ti){l-Tl) 

Note that D takes some finite value at q = 0, proportional to the residue of the pole from 
the massless graviton mode. Now we can compare the two quantities 



S)|g=0, 2)|g^-oo. (3.5) 

Suppose we find values for the input parameters such that there are no tachyons in the KK 
graviton spectrum. Then it must be the case that 

2)1^=0 X Slg^-oo > . (3.6) 

As we vary the input parameters, we may cross the hypersurface in the parameter space 
defined by £)|g=o x !l)|g^_oo = 0. Just across this boundary are models which each contain 
a single tachyon. To find this boundary between models with no tachyon and models with 
a single tachyon, we need explicit expressions for 3|g=o and £)|g^_oo- 
2)|g=o can be worked out straightforwardly: 

S|,=<. = -2E(^'°«5i^ + ^)=-2Cf , (3.7) 

where we used tanh~^ ^ = i log 

Before trying to evaluate 5)|g_+_oo) note that our solutions of the equations of motion 
should be real. While P'^^^^^^^^-^^^^i^z) is always real, Q'\^_y^^ ^j^^^ '^ j^k^") becomes a complex 
valued function for q < —9/4. Thus for q > —9/4 we still write ( p.llD in the form 

y^'^y) = ^^M+v/9T4^)/2(^) + ^Q(-l+^)/2(-) > (3-8) 

but for q < -9/4 we use relation 8.843 of and relations 3.6.1(4)-(5), 3.3.1(7) of |ll to 
write the real solution 

y'^Hy) = Apr\/.i^) + ^^^(/ + ^)(/ + l)Pi\M-^) > (3.9) 
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where we have defined 2ip = y/9 + 4g. The two expressions match at p = 0, q = —9/4. 

Plugging (^) into the brane-boundary equations of motion, we obtain ao, bo, ai, hh 
which are real for any large and negative q, i.e., large and positive p: 



bo = { - voip' + l)^/^iPr,\/2iTo) + 2^;\/2(ro)} , 

= { - vUp' + \)^^^lP-\/2iTL) + 2P7^\/,{Tl)} , (3.10) 

Using relations 3.2 (14), 3.4 (1), 2.3.2 (17) from [^], we can derive the asymptotic behavior 
of the conical function PipTi/2i^'^^^) large p: 

and thus we get 

2)1,^-00 « -4, ■ , (e^(^-^°"^^) - eP(^o+^-)) , (3.12) 

27rpv sm 6*0 sm 6*^ V / 

where cos 6i = Ti. 

Now we are ready to look at (|3.6| ), which implies 



2)|g=o>0 and 2)|g^_oo>0, (3.13) 



or 



2)|g=o<0 and 2)|g__oo<0. (3.14) 

Note that the sign of 3|g^-oo is the same as the sign of 

voVl{cos^^ To + cos^^ Tl — tt) . (3.15) 

Figure ^ shows the tachyon counting for models in the {vl,wl) plane when vo > 0- 
The sign of S)|g=o flips when we cross the line Cg^^ = 0, which is shown as the solid curve. 
The sign of ( |3.15| ) flips whenever we cross one of the two dashed straight lines. The vertical 
dashed line is just vl = 0. The slanted dashed line is 

WL = -6il-T^)vL-l2To . (3.16) 

Thus the shaded region in Figure ^ gives models which have a single tachyon (or are not 
AdS5/AdS4). The models in the four unshaded regions have an even number of tachyons. 
Models in the large unshaded region to the right of vl = are tachyon-free, as expected. 
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Figure 6: !D|g=o x £i|g^-oo > on (w^, u;i)-plane when vq ~ 2 and wq — I. The solid curve is 
where S|g=o = and dashed straight lines are X)|g^_oo = 0. 

By further numerical analysis, we find that models in the upper two remaining un- 
shaded regions are also tachyon-free, while models in the lower one have two tachyons. 
Once we also require Cg^^ > 0, only the middle unshaded area remains. This region is 
rather special. It contains no tachyons and obeys Cg^^ > 0. However in this region one can 
see numerically that the massive KK graviton modes are kinetic ghosts. 

Combining all of the above results with the additional requirement that > 0, i.e. 
no ghost radion, we get figures like Figure |^ for {vq,wo) chosen such that Tq > 0. Models 
which are ghost-free and tachyon-free correspond to the unshaded region in the {vl,wl) 
plane. For < 0, a similar analysis finds no region which is free of tachyons and ghosts. 

These results fit qualitatively with our physical intuition. Models with negative Mj 
are excluded. Models with too much negative tension branes are also excluded. 

A special case is that of an infinite extra dimension. In our general framework we 
make models of this type by sending L ^ oo. We do this by fixing wl = +12 and (for 
simplicity) vl = Q. In this limit the graviton zero mode is not normalizable and drops out 
of the theory, so there is no Cg'^^ to consider. As for C^, since Tl = +1, (|3.3|) becomes 



which is positive for any vq and wq < 12. Our tachyon analysis reduces to examining ao|g=o 




3x 



C^H^ {l + To){2 + vo + voTo) 
2k 1 + vqTo 



(3.17) 
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Figure 7: Allowed region (unshaded) on (w^, WL)-plane with vq = 2 and wq = I- 



and Oo 



— oo 



. Since 



ao|g=o 



(i-r2)<o 



and 



1 - eP^""-^"^ 



(3.18) 



— oo 



-Vo 



V27r sin 6*0 ^/p 



only t;o > is allowed. So in this class of models the allowed region is f o > and wq < 12. 
4. Strong coupling 

For generic choices of the four input parameters vq, wq, vl and wl we will not have a strong 
coupling problem, as long as we choose parameters away from the borderline between ghost 
and non-ghost regions, and away from limits where H ^ 0. 

To determine the strong coupling scale, we need to calculate the cubic action for the 
radion. See Appendix ^ for the details of the calculation. The general cubic expansion of 
the radion part of ( |2.1| ) is 





S 




2M3 ^3 
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(3) 

When it comes to the strong couphng, would be expected to give the most stringent 



limit. However with F{y) determined by ( |3.1| ), direct calculation shows that 



2 



+ TTr 1 , ^ g(^) + 



y=yi 



0, (4.2) 

,(3) 



The next candidate, Q , turns out to be 



'2 - /-'^^l 12 + fc2 1 + z 



gao-OL/ 2 , , 2 , -ISolTo"^ - 12aor^^ + 6(ao - aL)log4 



y=y^ k Jq 1 + Z 



= X — 7^ — |ao + «oaL + aL + i i 

where T^^ = tanh~^Tj and Li„ is a polylogarithm function. This does not vanish. Then 
from the full radion action 



(2M3C^)3/2 
we can determine the strong coupling scale: 



+ (^^3c W2 j d'x^gi^{{V^f - V^VW^V^V>} + • • • , (4.4) 



^1/2 

As an example, let's consider a DGP-like limit [v^ = 0, wl = +12, t;o is large). Then ( [4.5| ) 
becomes 

Asc ~ —^12-WQ, (4.6) 
Mo 

where our parameter Mq is equivalent to the parameter Mp in DGP. This agrees well with 

Hi- 

In our general framework we can understand the robustness of the strong coupling 
problem. As long as we restrict ourselves to models in the ghost +tachyon free region of 
the parameter space, the coefficient of the graviton kinetic term never vanishes, as can 
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be seen from ( 2.31 ). From ( [4 .51) we can see that the radion becomes strongly coupled in 
any limit where H ^ 0. This includes the DGP-hke hmit just mentioned, as weh as the 
"bigravity" hmit Tq +1. Unfortunately these are precisely the limits in which we find 
ultralight graviton modes. 

5. Green's function analysis 

In [^] , the set of coupled equations of motion of the graviton and the radion was obtained 
in a straight gauge. Once we eliminate the radion, we get four independent equations 
involving the graviton only. These are of the form: 

^iV^'"^ = 0' (5-1) 
Pg/i^'^ = 0, (5.2) 

■^fiupa"' \y=yi 



rtp'>n.=.. = 0, (5.3) 



where the full expressions are given in ( [B.6D -( |B^ . These imply the following Green's 
function equations in the straight gauge: 

7^(i) CP'' - ^ n . :(r y /^'H^ - a:0<5(j/ - 

^^^G|;V=0, (5.5) 
[^?SGf;,4=..=0, (5.6) 

where the AdS^ bitensor 0^y-^iyi{x^x') is given below. Then, for any given source, we get 
the linearized solution for the graviton from 

h^uix,y) = / d'^x'dyy/^Gf,^.^'u'ix,y;x',y')T'''''' {x',y') . (5.7) 



In Appendix the Euclidean versions of (p^)-([5^) are explicitly solved to obtain the 
Euclidean Green's function. After dropping 4d total derivatives which will vanish when 
contracted with a conserved stress tensor, we can write the following expression for the 
Euclidean Green's function on the 0-brane: 



Gf,u;p'u'{x,x',y = y' = 0) = Y^Gi{u;pj)l T^J.^,^, 

j 



where the sum is over the residues of poles from the individual KK graviton modes. The 
masses mj of these modes are determined by the solutions of ( p.l2| ) , and are expressed in 
terms of the parameter pj: 

m] = H\p] - h . (5.9) 
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The variable u is related to the geodesic distance // between the points x and x' in AdS4^: 
u = coshff/i — 1. The function Gi{u;pj) is given by 



1 d 



Gi{u;pj) = — — 7-T-^(3p,-i/2(l + u) , (5.10) 



where Qp--i/2 is a Legendre function. The AdSi^ tensor structure of the Green's function 
( |5.8| ) is contained in the bitensors T^^^ and T'^'^\ which are part of the complete bitensor 
basis given in ( [8.121) : 

'^!n};,i'u' = d^d^/ud^d^'U + d^d^'ud^d^'U . (5.11) 



The bitensor on the right hand side of ( ^.4[) is given by: 



0^u;,i'u'{x,x') = ^{gfj,f,'gpu' + 9iiu'9uf,') , (5.12) 
where we have defined bivectors: 

9(1(1' {x,x ) = d^df^'U ■ (^-13) 

The bitensor 0^j^y^^iyi{x^x') is determined by requiring that it has the property: 

0^,,,^>,>{x,x' = x)T^'''''{x) = T^,{x) , (5.14) 
which can be checked easily using ( 5.12| ), ( |5.13| ). 



Apart from the j-dependent coupling constant m'^\ the expressions ( ^^ ) and ( 5.1C| ) 



are the same as those derived by Naqvi [ 43 1 , p^] for a massive symmetric tensor in AdS^. 
Note this is very different from the Green's function for the DGP model, which in addition 
to an overall coupling constant has extra gauge-dependent tensor structures |4^] which 
diverge as Mp oo. In our straight gauge analysis of warped DGP-like limits, such 
effects are completely absent. The breakdown of the linearized gravity approximation is 

due entirely to strong coupling of the radion, as we discussed in the previous section. 

(i) 

The coupling constant is given by 

1 _ bLP^'^^^,{-To) - aLQ% _^/,{-To) 

. (5.15) 



Mp)2 2feM3 [9,(ao6L-6oaL)],=,^=p2_9 

This is the effective 4d gravitational coupling constant of the j-th KK mode of graviton to 
matter on the 0-brane. Depending on the choice of input parameters, this coupling may 
show crossover behaviour; if the values of ( |5.15| ) for modes heavier than some mass Aco are 
highly suppressed compared to lower lying modes, then Aco defines a crossover scale. 
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5.1 asymptotic behavior of the Green's function 

The u ^ oo limit of the graviton Green's function ( |5.8| ) shows how sources on the 0-brane 
interact for geodesic separations which are large compared to 1/H, the AdS4 radius of 
curvature. The asymptotic formulae can be extracted from [41|, giving: 

= y' = 0)^^^ Y.Gii^'PM ^S.v - r^^^'T^'J;,'u'} > (5-16) 

where 

4,2^2^0)2 rfe + i) 

For large u, u ^ ex.p{Hfi). Thus from ( 5.17| ) we see that the contribution of each graviton 
mode to the Euclidean Green's function damps exponentially for large geodesic distances. 
In models with a g < —9/4 tachyon, pj is imaginary, and the asymptotic contribution to 
the Green's function oscillates, as for a flat space tachyon. 

5.2 vDVZ discontinuity 



To investigate the vDVZ discontinuity problem, we will follow j46]-[^]. In the flat 4d 
spacetime limit {Hjj, ^ 1, i.e., u — > 0), expanding the Euclidean Green's function 
gives 

G^y-^'y'{x,x',y = y' = 0) 

+ ^f^u'^u'u' - ^ • ^^-^^/luSf^'u'^ ■ (5.18) 



As qj varies from to 00, ( 5.18| ) varies smoothly between the flat space Euclidean tensor 
structure of a massless 4d graviton: 



CP^^'v fx —{^^I.^l'^vu' + S^f'Su^i' - Sfj.uS^'u') , (5.19) 



and that of a massive 4d graviton: 

C^if-^V^ Of —i^fifi'^uu' + S^lu'SufM' - ■::Sf,uSfi'u') , (5.20) 

The condition for an ultralight graviton mode to avoid the vDVZ discontinuity problem is 
that Qj — > in the limit — > 0. 



6. Warped models with an infinite extra dimension 

As an application of the different results developed in the previous sections, we will highlight 
here the phenomenology of a model with an infinite extra dimension. This generalizes the 
single brane Karch-Randall (KR) model to include a localized curvature term. In order 
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to get L = cxD we fix wl = 12 (thus = 1). Also we choose for simphcity vl = 0. The 
model is therefore completely determined by the higher dimensional Planck mass, M, the 
bulk curvature k, and the two brane parameters vq and wq. The ratio of brane to bulk 
curvature reads 

H 



X = Vl-^o, (6.1) 

where Tq is given in terms of vq and wq by equation (^.6|). Different limits of this model 
have been studied by a number of authors |T9|],|49|-H|. 



The infinite extra dimension results in the graviton zero mode being non-normalizable 
and therefore it decouples from the spectrum, whereas for strong warping (i.e. large k/H) 
there appears an ultralight (mi ^ H) graviton mode. The ultralight mode couples to 
brane matter with comparable strength as the RS zero mode, while the rest of the graviton 
spectrum is heavier than H and has much weaker couplings due to the warping. For large 
localized curvature term, the mass of the different modes decrease as we increase vq, with 
the first mode again becoming ultralight. The coupling of the KK modes to brane matter, 
on the other hand, is suppressed only for modes heavier than some crossover scale. This 
crossover scale depends on the size of vq. In the following we shall see how our general 
results, when particularized to this model, reproduce these features. 

All our general equations can be readily adapted to an infinite dimension by noticing 
that Tl = 1 implies = 0. The y-dependent part of the graviton KK mode wave function 
reads now 

^^^^(-) = T-i+v^)/2(-)' (6-2) 

with z = tanh A;(y — yo) and ky^ = tanh""*^ Tq. The masses are given by rin? = qH^, where 
q is determined by the zeroes of the equation 

ao{q) = 0, (6.3) 



with oo given in ( |2.13 ). The kinetic coefficients for graviton KK modes and the radion 



can be obtained directly from ( 2.32| ) and ( |3.3D by taking the corresponding values of the 
different parameters, 

= ^{'^ j\ y^''^^dz + VQ{l-Ti)y^'^'^\-TQ)) , (6.4) 

_ 3x'gMl + ?o)(2 + ^o + ^oTo) 
^ 2k 1 + vqTo ■ ^ ' ' 

We have not written cf'^ because, as expected, Tl = 1 renders the massless graviton 
non-normalizable and therefore it decouples from the rest of the spectrum. Similarly, the 



couplings of the KK modes to the brane at y = can be obtained from ( 2.36| ). As 
an example, we show in Figure ^ the masses of the first five modes (in units of the brane 
curvature H) as a function oi k/H, for vq = (this is the original KR model). The warping 
makes the mass of the different modes smaller as k/H increases, and in particular singles 
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Figure 8: Masses (in units of the brane curvature H) for the first five massive graviton KK modes 
as a function of the bulk to brane curvature ratio, k/H. The strong warping regime corresponds to 
the far right region of the figure. 
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Figure 9: Masses (in units of the brane curvature H) for the first five massive graviton KK modes 
as a function of vq for fixed H/k = 0.25 (left panel) and H/k = 0.95 (right panel). 

out the first massive mode as an ultralight one, with mass much smaller than the brane 
curvature. In Figure ^ we show the effect of the localized curvature term by plotting the 
masses as a function of vq for fixed value of H/k = 0.25 (left panel) and H/k = 0.95 
(right panel). The effect is clearly not as dramatic as with the warping as we can see by 
comparing the slopes of the curves with the ones of Figure ^ or the one on the right panel 
(small warping) with the one on the left (large warping). 

Now consider the couplings of the graviton modes to brane matter. In Figure |l^ we 
show the couplings of the different KK modes to the brane at y = in units of the first 
mode coupling. Again the warping effect is clear, with very suppressed couplings for large 
warping (large k/H). In Figure |ll|, we show the coupling as a function of fo, for fixed 
values of H/k = 0.25 (left) and H/k = 0.95 (right). In the left panel we see how the 
warping makes the effect of the curvature term less acute, whereas the right one, where 
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the warping is very small, clearly shows the effect of vq. This plot also shows that the 
couplings of the different modes are more or less suppressed depending on their masses (a 
clear hint of the crossover behaviour). 
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Figure 10: Coupling of the first five massive graviton KK modes to the brane at i/ = (in units of 
the coupling of the first mode) as a function of the brane to bulk curvature ratio, H/k. The strong 
warping regime corresponds to the left of the figure. 



Figure 11: Coupling of the j — 2,3,4,5 massive graviton KK modes to the brane at y = (in 
units of the coupling of the first mode) as a function of vq for fixed H/k = 0.25 (left panel) and 
H/k = 0.95 (right panel). 



The crossover behaviour is best seen by plotting the coupling of the different modes 
as a function of their masses, for different values of vq. We do this in Figure where 
we have fixed H/k = 1 (small warping) and the different sets of points correspond to 
Vq = 0.05, O.f, 0.15, 0.25 and 0.4, from top to bottom. We have superimposed (solid lines) 
the expectation for these couplings in the DGP model as given in [pO| , with a crossover 
scale Tc = vo/k. It is clear from the figure that once the localized curvature term is large 
enough to suppress warping effects the agreement is quite good. 
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Figure 12: The dots represent the coupUngs of the massive graviton KK modes to the brane at 
?/ = (in units of the coupling of the first mode) as a function of their mass (in units of H), for 
fo = 0.05,0.1,0.15,0.25 and 0.4, respectively and fixed H/k — 1. The lines correspond to the 
expectation from the DGP model with crossover scale rc — vo/k. 



Let us now turn to the Green's function analysis in the KR model. Using the fact that 
cl = we can write (|B.67 ) as, 

Gi{u,y = y' = 0) = 



dp ■ 



-1 



Then we immediately see that our spectrum is determined by solving 



(6.6) 



0, 



which of course gives the same result as the KK analysis. We now contour-integrate (| 
to get 



(6.7) 



0) 



Qi{l + u) 



Gi{u,y = y 
1 

where the sum is over the poles of ( |6.7D . Comparing this with ( |5.1C1| ), we obtain 

1 H r {v+\)p-\-n- 



2M3 



(6.8) 



(6.9) 



As we show in Figure the couplings computed this way (dots at y = in the figure) 
agree with the ones computed in the KK analysis (lines for couplings of the j— th KK 
modes, with j = 1,2,4,6 to a probe brane as a function of its location). We have fixed 
H/k = 0.95 and f o = 1 in this plot and consider values of y < 2yo- It is also trivial to 
check that the masses and couplings agree with the results in |3S], by taking vq = in the 
equations above. 
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Figure 13: Coupling of the j— th graviton KK mode (in units of H/M^) to a probe brane localized 
at ?/ for < ?/ < 2?/o- We have fixed H/k = 0.95 and = 1 and j — 1, 2, 4, 6 from top to bottom in 
the coupling at the brane at y = 0. The dots correspond to the couplings to the brane as computed 
in the Green's function analysis. 



7. Conclusion 



One of the frustrations of trying to understand braneworld gravity is the large number of 
models which display various combinations of interesting or problematic features. In this 
paper, following [^5|, we have addressed this difficulty by introducing a general unified 
framework. This framework describes general 5d warped setups with two branes and local- 
ized curvature. Apart from the 5d Planck mass M and the bulk curvature k, there are only 
four input parameters: {vq^wq^vliWl)-, related to the two brane tensions and the strength 
of the localized curvature. This framework includes as special cases both Randall-Sundrum 
models, both Karch-Randall models, the Lykken-Randall model, bigravity models, as well 
as a large class of DGP-like models, where the localized curvature competes with or dom- 
inates the warping. 

Our unified framework has the additional advantage that we can gauge-fix the 5d 
general covariance explicitly, and exhibit the gauge-fixed 4d effective action for the physical 
degrees of freedom. This we have done in section 2, building on the results of |35|. In section 
3, we used this effective action, together with an analysis of the KK graviton spectrum, to 
reveal the presence of kinetic ghosts and tachyons in a large class of models. Depending 
on choices of our four input parameters, we found examples of models where the radion is 
a ghost, the massless graviton is a ghost, or the massive graviton modes are ghosts. We 
also found models containing one or two tachyons in the graviton spectrum. 

Our framework is especially well-suited for examining the problem of strong coupling. 
In the literature, the phenomenon we call strong coupling has shown up in many confusing 
guises, all having to do with a breakdown of the linearized approximation for braneworld 
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gravity. Here, we mapped out the parameter space in which either the radion or the 
massless graviton become strongly coupled. We found that the graviton only becomes 
strongly coupled for models which also contain tachyons and ghosts. We also noted that, 
in our straight gauge formalism, the graviton propagator does not have tensor structures 



which diverge as an ultralight mode becomes massless [45]. Because we are always in a 
straight gauge, we also don't have to worry about large brane-bending effects ]55[ ]. 

Thus in our framework strong coupling is an isolated property of the radion, and can 
be studied in a straightforward way. In section 4 we have written a general formula for 
the resulting strong coupling scale in our framework. We observed that there is a tight 
relationship between strong coupling and the presence of an ultralight graviton mode. This 
was known for the DGP and bigravity models, but here we see it in a unified picture that 
interpolates smoothly between both kinds of models. As expected, the vDVZ discontinuity 
is never present for ultralight modes in our warped setups. 

In section 6 we studied the effects of localized curvature in the subset of models that 
have an infinite extra dimension. These are the models in which the usual 4d massless 
graviton is absent, having decoupled in the limit that we removed the second brane. We saw 
that localized curvature and warping both have the effect of making an ultralight graviton 
from the first massive KK mode. We compared these two effects in models which contain 
both. The dimensionless parameter vq controls the strength of the localized curvature, 
while the dimensionless ratio k/H measures the strength of the warping. By looking at 
models with comparable values for these parameters, we saw that the warping effect is 
stronger. Thus "locally localized" gravity is more efficient than localized curvature in 
creating a simulacrum of 4d gravity. 

Finally, in section 6 we exhibited explicit crossover behavior for models in our general 
setup. In these models the couplings of the KK graviton modes to brane matter are 
unsuppressed up to a mass scale l/vc ~ k/vo, and are highly suppressed for modes heavier 
than this scale. This behavior is the most interesting phenomenological feature of DGP 
braneworld gravity; here we have found it in a general warped framework. Our results, as 
seen in Figure |l2|, agree quantitatively with what one would predict by analogy with DGP. 

We have opened a window to a new class of models which have crossover scales from 
four dimensional gravity at short distance to five dimensional gravity at longer distances. 
We have presented a well-defined framework to analyze strong coupling behavior, which 
remains a fundamental obstacle towards developing realistic effective theories of braneworld 
gravity. 
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A. Radion action to cubic order 



With /i^4 = 0, the general cubic expansion of (2.1) is 



S 



2M3 



cubic 



3a'2 3 -3if2 + i2a^ 9a' r/, 2 
t^iIaa H r — 5 ririAA H — r — il iIaa 



oCL 8 o 4(X 



-^/l44V^(/i'^"V,/l^) - ^/l44V,/l^"V^/l^ - ^/l44/i'^°V,V^/l^ 

.J a a a 

^ ^ -I y=yi 



(A.l) 
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The following identities will be useful: 

a ^ a 

^ Qj Zi Qj 

1 2a' 

Upon making the substitutions 

V ^ a^3^i(y)V^V^i/; + gi,uy2{y)4^ , = , (A.2) 



( A.l ) produces a few hundred terms. Collecting terms with three 3^i's, the bulk part gives 



it2 q ij2 o i:j'2 

Performing 4d integration by parts a few dozen times, a dramatic simplification occurs and 



|A.3|) becomes: 



■ (V'(V V)' - V'V^V^^^V^V^V - ^V^'V V) } • (A.4) 



f2 

i2 



By a similar calculation, we can show that no term with three 3^i's survives in the brane- 
boundary part. Note that terms with eight V's are all cancelled. 
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Repeating the above procedure with the remaining terms, we get the radion cubic 
action: 



S 



2M2 



3 9ij2 

+_(__V,(vW + ^*^v=*)} 



+- 



16 
. 3 



16 



o 



6i?2 X „/6a'. 



2p2( 



3a'2 9a' ,\ 



-3i/2 + i2a'2 



3a' 



/2 



+ 



+- 

/3if 



3 ,,-2 m2 9//^ 



Collecting terms of the same V-structure, we get ( [4.1D . 



(A.5) 



B. Graviton Green's function 



In this section, we explicitly calculate the graviton Green's function in the Euclidean sig 
nature. 
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In |35], the bulk EOM for the graviton and the radion are obtained to be 








3H^ ~ 

2 dfivh 



-n-fj^u + 9^lu^^ + ~9fj.un H rifiu 

3a' 



-VuVuh^"" + V^/i + — h 
a 



3a' ~, m'^ 



h 



12a' 



12 



a' a' 
(a^h!)' ~n , 4a' ^, , , , ^ , 



a 



(B.l) 
(B.2) 
(B.3) 
(B.4) 



and the brane-boundary equation is 



i{hfj.u' - g^ivh') + 



(B.5) 



Using (|B.2| ) or ( [B.3| ) we ehminate the radion and get the equations involving the graviton 
only. The bulk part gives only two independent equations: 



,// 1?// a 7/ 4a'% ~ ,^ ^ 

-f^t^u + 9iiv-:h + —9fj,iyh H g""/^^ + T^9fiiyh (B.6) 

+^(^(v^Vp/ie' + V.VpO - V^V./i') + 5^,^(v2/j' - VpV^hn , 



SH"^ ~ ~ 4a' ~ 



i2 M 



(B.7) 



and the brane-boundary part gives 



+ 2kTijhf,^-^(Xi + 



(B. 



Then, from 



and (B.7) we get the equations which the graviton Green's function 
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should satisfy, whereas ( [B.8D provides the boundary conditions on the branes: 



2 ~ ~ „ 1 ~ 2 

1 a' 4a' ^ i?^ 



(B.9) 



1 . 



2M3 



^9 



5{y-y'), 



a ~ 



(B.IO) 



+ 2A:rj ) G^y-^iyi - ( Ai + -^^jj^ ) g^iuG'^p.^,^, 



2/cr,- 



A. 



1 



(B.ll) 



Now we fohow Ql and |Q to solve ( B.9p -( prTl ). First we decompose Gf^y-^j,iyi using the 
five independent bitensor bases: 



G,u,^u' = Y.T'^^p'u'A^^^y^y') 



i=l 



where u = cosh. Hfi — 1 with /x the geodesic distance between x and x', and 



rp(2) 

7^(3) 
7^(5) 



9pv 9p'v' ) 

d^uduud^'udiy'U , 
d^d^iudydyiu + d^d^mdyd^m , 
g^ud^'udy'U + Qf^iyid^uduU , 
d^dfj^iudyudyiu + dyd^iud^udyiu + (/i' 



Then we reorganize (B.12) into 

G^.;pV = r(3)a4Gi(n,y,y') + r^'^G2(u,y,y') 



(B.12) 



(B.13) 



(B.14) 
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where 

K;,i'v' = gfj.'u'duUa'^A{u,y,y') + d^>udu'uduUa'^C{u,y,y') 

+du{d^>udy>u) a^B{u, y, y') , (B.15) 

^^lI^■,u' = OiMud^'U a^A{u, y, y') + d^uduud^iua^C{u, y, y') 

+d„'{dfj,udiyu)a'^B{u,y,y') . (B.16) 

With this reorganization, we can single out the physically meaningful part, Gi and G2. 
The ^W's and Gi,G2,A,B,C are related by 

A« = G2 + 4i72(i + n)A, ^(2)^4^4^^ A^^^ = a\Gi + 4B) , 

A'^^^ =a'^f = a'^{2A + 2H^{2B + {l + u)C)}, A^^^ = 2a^ {B + C) , (B.17) 

with an overdot implying d/du. Upon plugging ( |B.14[ ) into ( |B.9| ), we get 



1 (^(4) {^j. _ 2;') 

-Xa^l^l'9vu' + 9f^u'gufi') p - v') 



2M^H^\{2 + uf 2 + V a4 47r2M 



(B.18) 



where we have converted + Ofj^u'dufi' into T^^^'s applying the relations given in |44|. 

To convert 5^^\x — x') into S{u), we start from a Euclideanized global coordinate system 

(t, p,0,4>) 

-^ds"^ = cosh^ p dr"^ + dp^ + sinh^ pldO"^ + sin^ 9 dcf?) , (B.19) 

where u measured from (0, 0, 0, 0) to (r, p, 0, 0) can be explicitly written as 

n = cosh/ocoshr — 1 . (B.20) 

Note that at a given n, r is restricted by — cosh~"'^(l+u) < r < cosh~^(l+ti). Transforming 
to (t,u,6',0), 

, n ii(l + 'u)2(2 + 'u)sech2r , n 1 ,n 
^ds^ = '-^ '-^ dr^ + ^du^ 



a 



^ (1 + - cosh^ T (1 + n)2 - cosher 

2(l+u)tanhT , , (1 + n)^ — cosh^ r , o „ , ,0. ,^ 

^ ' dTdu + ^ ^ — 2 {dd^ + W OdctP') . (B.21) 



(1 + w)^ — cosh^ r cosher 
Then, when integrated with a function depending on u only, 
5^^\x-x') 6{T)6{u)6{e)6{(l)) 



^ ^ (1 + u)sech?T\J (1 + n)2sech^r — 1 sin 9 

6{u)6{t) 



" 47r(l + 'u)sech2T y (1 + 'u)2sechV - 1 
ff^ 6{u) 



(B.22) 
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In the last step we use, for arbitrary function / of u, 



[ dudry'^i^^P^fiu) = /(O) = / dudr^^^j^fiu) , (B.23) 
J yJg{T, u) J J dT^ 



and then replace 5{t)/ y/g by 1// dr^/g, where 

<i'T\/g = / dr (1 + ti)^sech^TW 1 

7-cosh-i{l+«) V (1 + ^)' 



tanh^ r 



|u(2 + u). (B.24) 



p3q) gives 

g^'u'duU -E5 + dfj^iudyiudyU Eq + dy{d^iudyiu) Ej = . (B.25) 
Using several identities presented in |44],the Ei^s can be worked out: 



4 



+^^(2 + 'u)G2 + — 7^i(2 + u)G2' 

+^(1 + n)V + ^u{l + n)(2 + «)/ 

/ 3iJ2 6i?2 4ij2 3q/2 2ij2 . 
+ ^u{2 + u) ^ 1 + uY + + 2 + n / 1 + u 2 + n / 

f a H'^ 2a' , , a -ff ^ \ 1 q 772 

/24if2 16^2 ^ 2ai/2 

+ — ^ 1 + -u r + — 2 + n M + — ^-u 2 + u)A' 

, ,. a if2 



-n 



(1 + n)(2 + + -— «(! + 'u)(2 + u)A' 



a2 V a a' a''' J 



-^(1 + u)C-^^{l + u)C', (B.26) 



ii/o — Z H^G\ — — -H^Gi 

3 3a' 



80, , ••• 8 n ■• 4a n , , ■■■/ 4a n ■•, 

--H'il + u)B +-H^B - —H\l + u)B +^H^B' 

+ ^H\l + U)C + _ 48„/2^J^ _ ^^2^// 



-33- 



E3 = Z+ {4H^ - 12o'^)Gi - ^H^{1 + u)Gi + (- 8aa' + ^H^)Gi - H^u{2 + u)Gi 



-48a' 



a 

3^ 

'^B + ^^H\l + u)B + UM'^ + u)B + (^4H' - ^a') B' 



-4a^B" + ^H\l + u)B' + ^hM"^ + u)B' 
da' 6a' 

+-H\l + u)C - -H\(2 + u)C + ^H'^il + u)C' - ^H\(2 + u)C', (B.28) 
3 3 3a 3a 

E, = {-2H' + 6a'') Gi - 2H\1 + u)G, + (—a' - ^H^)g^ - \h^G, + \a^G^" 
^ ' \ a 2a! J 6 2 

2a! ^ ' 6a' 
+10(l + u)/ + t/(2 + u)/ + 8/ 

+10/ + + 14(1 + u)i + ^(1 + n)f + 2u{2 + u)] + ^u{2 + u)f 



+ (24a'' - 72ij2)5 + (l^a^ - + 2a^B" + ( - 104/f2 + 24a") (1 + u)B 

- (-H\l + uf + ^HM'^ + u))b+(- + ^a^) (1 + u)B' + 2a\l + u)B" 

16a' o 4a 



+ ( _ 88i/2 + 24a") (1 + «)C + (^a^ - ^/f^^ (1 + u)C" + 2a2(l + 
+ ( - ^H^u{2 + -a) + 12a'\(2 + n) - ^i?'(l + u)') C 

\ 3a' 3a' a / 

-8i7\(l + 'a)(2 + u)C - ^H^u{l + n)(2 + u)C' , (B.29) 



E, = -^Gi - ^H'G,' - 4G2 - ^G2 + ^H'f + 2H\1 + u)f 
3 3a' a' 

8a' 



-I2a''f+2H\1 + u)f+ (-^a' + ^hA/' + ^H'{1 + n)f - a'f" 

\ a a' y a' 



-8H'^{{1 + u)A + 2A) - ^H\{1 + u)A' + 2i') 

a' 



.^B - P-H'& + + P-H^C' , (B.30) 

3 3a' 3 3a' ' ^ ' 
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i 6 6a' ba' 

+ C4h' - 24a") B + IhHi + u)B + |i?\(2 + u)B +C4^H^ - ^a^)^' 
V 3 / o 3 V 6a' a J 

-2a^B" + p-H\l + u)B' + ^H''u{2 + u) B 
6a' 6a' 

+ (^H' - 2Aa")c+ -H\l + u)C - -HM^ + n)C + ( - ^a')c' 
V3 / 3 3 VSa'a/ 

-2a^C" + ^H\l + u)C'- ^H^u{2 + u)C', (B.31) 

Eg = -(6ii"^ + AH'^a'^)Gi + lOH'^il + u)Gi + 2H\{2 + u) Gi -2aa'H'^Gi 
+18(1 + u)G2 + 3«(2 + u) G2 

-36(1 + u){H^ - a'^)f -3{H^{A+18u + 9^^^) - 2u{2 + u)a'^)f 
-3ij\(l + n)(2 + u)f+ ISaa'il + u)f' + 3aa''u(2 + u)f 
+144i?2(l + u)A + 36H^{2 + 6u + 3u'^)A + 12H^u{l + it)(2 + u) A 
+16H^ ( - 3^2 + a'^)B + Saa'H^B' 

+ ( - 24i7^ + 32H^a'^)C - 2AH'^{1 + u)^ + leaa'F^c" ^ (b.32) 

42H^ - Aa^)Gi - 6ff^(l + u)Gi - 2H^ Gi -2adGx 
+42 f + 30(1 + n)/ + 3ti(2 + n) / 

+ ( - 336i^2 + 160a'2)^ - {240H'^ - 16a''^){l + u)B - 24H^u{2 + u) B 
+80aa'B' + 8aa'{l + u)B' 

+ ( - 168H^ + 176a'2) C - 8(1 + fx) {blH^ - 28a'^) C 
-12(if2(l0 + 26u + ISn^) - 2u(2 + u)a'^)C - 12H^u{l + u){2 + u)C 
+88aa'C' + 112aa'(l + u)C' + 12aa'u{2 + u)C' , (B.33) 
£^8 = + a'^)Gi + 2{l+u){-3H^ + 4a'^)Gi 

-2H^Gi + 18aa'Gi + 4aa'(l + u)Gi 
+18/ + 24(1 + u)f + 3n(2 + u)f 

-U4{H^ - {a'f)B - 16(1 + u){12H'^ - 7a'^)B - 8u{2 + u){3H^ - 2a'^)B 
+72aa'B' + 56aa'(l + u)B' + 8aa'n(2 + u)B' 

-8(1 + u){21H^ - 10a'2)C - 8(3^^(4 + lO^x + 5u^) - u{2 + u)a^)C 
-12H^u{l + u){2 + u)C + 40aa'(l + u)C" + 4aa''u(2 + u)C' , (B.34) 

with 

f = -2A + 2H^{2B + {l + u)C), (B.35) 

Z = 2f - 2f - ^f' = -8A - ^{A' + H'^{2B' + {I + u)C ')) . (B.36) 
a a 
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The brane-boundary part can be worked out similarly: 

+gt.ud^,^ud,,u H'^Ef^'^ + g^.,.d^ud,u + T^^'^ Ef\=y^ = , (B.37) 

where 

. _ {^n^,T. -^ ^ ( A ^ 4A.) ) + ^ + 2A.) (1 + u)G. 

+^(1 + u)V + ^u{l + u){2 + u)f 

- {^mu{2 + u) + ^{l + 6u + 3u^) ( + A,) ) / (B.38) 

+^ + 2Ai) (2 + 6tx + 3tx2) i + ^ + 2Ai) n(l + w) (2 + 

-4il^ (^T. + ^ (4 + 6A.) ) 5 - 2e.H^B' - ^ + 2A.) (1 + u)C , 

^) = -XiH^Gi + Xj+SkTiG^C + 4^ia2(7' , (B.39) 
^(bb) = ^( _ ^^2^(2 + - 2H^{1 + + AH^Gi) + 2kTia^Gi + ^ia^Gi' + A^/ 

+8kTia^B + 40^0^5' , (B.40) 
£;r^ = - (lOA,ij2 + ^ + a^A^T,) Gi - (2A,iJ^ + (i + 

V2A;rj 7 ' 2 ^ 
+4Ai/ + 5Ai(l + u)f + ^u{2 + tx)/ 

+K 2^ + ^0 ^" + ^< 2^ + ^0 + ")/^ + K 2^ + + 

- (72^2 + A,) + 4a''kT^ B - 2eiC?B' - (96H^ + A^) + Aa^kT^'^ (1 + u)B 

-2eia^{l + u)B' - 12H^ + \i^u{2 + u)B 

- (SAH^ + Ai) + 4a2A;ri) (1 + u)C - 2dia^{\ + ^x)C' 

- (\2H^ + Ai) (m(2 + + 4(1 + + 2a'^kTiu{2 + «)) C 

-6ij2 + Ai) u(l + n) (2 + -u)^ - 0ia\(2 + -u)C" , (B.41) 
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^4^^'^) = -XiH^Gi - XiG2 + 2kTia^f + eia^f' + 2XiH^f + XiH\l + u)f , (B.42) 
4^^) = 4kTia\B + C) + 2e^a\B' + C') + ^{H\1 + n)Gi + 2H^Gi) + A,/ . (B.43) 

£^2 5 and -E5 give the bulk equation for Gi: when u or y ^ y' , 



/ du' / dn"£;2 + 2 / du'E^ - E3 
J 00 J 00 J 00 



H\2Gi + 4(1 + u)Gi + u{2 + + 12a'^Gi + Saa'Gi ' + a^Gi " = , (B.44) 



and for li = and y = y' , 

E2 + 2E5 -E3 = (pi44|) 
1 1 



87r2M3a2U2u li- u-' 



6iu) + ( - - + ^)<5(n) - -5iu)}5iy - y') . (B.45) 



Similarly from £^2'^'^'*, E^^^ and £5'^'^^ we get 



du I du"E^^^^ +2 1 du'E,^ 
X H'^ 

- 2kTia^Gi + ^-—(4(1 + u)Gi + u{2 + u)Gi) - OiU^Gi ' 



,,p(bb) _^(bb) 



0. (B.46) 



y=yi 



Note that when performing integrations over n, we set the integration constants to be zero, 
because we want our solutions to die off as u gets large. Assuming 



Gi{u,y,y') = j dp 



X{x;p) 



with X = 1 + u and z = tanh — yo), ( |B.44 ) is written 



{l-z'yY{z,z';p)^ 



(B.47) 



d^x 



dX 



X 



= .(1-.Y{(1 



22;- 

dx'^ dx 
^ dz^ 



z^?Y 



dY 4 
dz 1 — z"^ 



2z^ - r—^Y^ =0. 



The complete orthonormal basis for X is given [56| by 



X{x;p) 



r(ip- 1/2)^1 



whose eigenvalue is + 4, i.e., 



d 



d 



(B.48) 



(B.49) 



i^-^):noX{x;p)-2x-Xix;p) 



p) = + 1) X{x; p) . (B.50) 



dx"^ ^ '-^^ fj^rf. 

Note that these bases satisfy the boundary condition at u ^ cxd because Plp_i/2(.^ + ^ 
£qj, ^g^j.gg ^_ Then for any given p, the solution for Y is 



ip-l/2V 



(B.51) 
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If we put the reference point y' between y = and y = L, we actually have three copies of 

dmi): 



(B.52) 



where Xi = (0, y'), X2 = {u' ■, L) and X3 = {—L, 0). Also the correct interpretation of ( p.4(^ ) 



IS 



I _ ^ + 9) (1 _ Tf) + 2m}Y{yi) - \{etY'{yt) + e-Y'{y-)} = , (B.53) 

where Of = —1, 9~ = \ ^ . Continuity of y at y = 0, L and y' gives three equations for 
AiS and Bis: 



(B.54) 



(B.53) gives two more equations: 



>1'(0+) - Y^'{{)-) + {kv^q{l - T^) + 4A;ro} ^(0) = 
Y^'{-L+)-Y^{L-) + {kvLq{l-Tl)+AkTL}Y2{L) = Q, 



(B.55) 



with q = —p^ — 9/4. The integration of ( [B.45| ) over {y' ,y'^) provides the last equation: 
X 



dp dl 

^Y^n(2 + n) 

4 



■Y' 



y=y^ 
y=y'' 







" ^^ r(ip-l/2)^;-i/2(l + ^) _ ^, y=y" 



T{ip) {n(2 + n)}3/2 



cosh fcun f/l 1 2\^,, / 1 2\-,, 



Using, 



r(ip + l/2-/i) 



r(ip) 



(B.57) 



and 



+ .1. .{J -] , m = 0,l,2,. 



m!r(z^ — m + 1) \2 



(B.58) 



^In jssj Si was introduced to represent the "outward" direction on the branes. That is, when we consider 
< y < L interval, on the 0-brane "outward" is in the direction of decreasing y {6q = — 1, where + implies 
we are on the right side of the brane), whereas on the L-brane we will leave the interval by moving in the 
direction of increasing y{9^ — +1). Then, if we are in the — L < y < interval, by the same reason we 
define 60 = +1 and 61+^ = 6+ = -1. 
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( |B.56|) becomes 

dpj^ dJ-^^^^^^P%, + n)4_i/2(l + u).r 

'\Yi{y'+-p')-Yi{y'--p')] 



y=y'^ 



y=y' 



T{ip' - 1/2) 


T{ip') 


T{ip') 


T{ip' - 5/2) 



cosh kyQ 



1 1 2 



+52(-iK2 + n)}3/2p3.^,_^/2(l + , 



cosh'^/cyo/ r(-ip' + 7/2)\ 
87r2M3 V r(-ip' -5/2); 



(B.59) 



or 



v'( '+ '\ vU I- i\ cosh ky^ i V{-ipl + 3/2) 

The general solutions for A's and S's with arbitrary y and y' are too lengthy to be written 
down. But since we are mainly interested in the gravity on the branes, we set y = y' = 0, 
to get 

where 



r(-ip) T{ip) aobi - boai 

Oj and bi are given in (p.l3|) with q and ^9 + Aq replaced by —p^ — 9/4 and 2ip respectively. 
We can decompose ( [B.61|) into modes by finding the poles of ( [B.62D , i.e., 

ao^L - boai = . (B.63) 

Once we choose our parameters from the allowed region, all the zeroes of ( |B.65 ) occur at 
non- negative q. 

Let's perform contour integration to evaluate the integral in ( [B.61| ). Using [41| 3.3.1(8), 

4-1/2(1 + U)='- COthTTp {Qjp_i/2(1 + U)- Qiip_i/2(1 + u)} , (B.64) 



( B.61 ) becomes 

Gi(^,y = ,' = 0) = ^^^-y^ dpp 

(B.65) 



bLP-\,^{-T,) - aLQ%.,„{-To) 
aobh - bottL 
Q!p-i/2(l + ^)-QVi/2(l+^) 
A/n(2 + u) 
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With H 8.737.4, 



(B.66) 



we can snow — - — aobL-boaL — even m p. inen 



Gi{u,y = y' = 0) 



1 if r . ^LPrl,J-To)-aLQl_,,J-To)Ql^,,Jl + u) 



dpp 



u 



d{-p) i-p) , ,^ . \ 

dpp -^=^= — . (B.67j 



Sir'^kM^ vr aofti - ^qOl \/u{2 + u) 

Since Qlp_i/2^^ + w) ~ e~^^^^^ for large p, we close the contour below. Noting that poles 
occur at positive g's, i.e., at pure imaginary p's, we finally get 

Gi{u,y = y' = 0) = '^Gi{u,y = y' = 0;pj) 
j 

dp >>i.Pi\,2(-^«^ - •^Mlp-iM-To) <3i_ 1/2(1 + 



Air^kM^ L dq dg{aobL - fcoai) y/u{2 + 



Next, £"2 and give 

duE^ — / du' / du" E2 

50 ./ 00 J 00 

= -H'^u{2 + u){B-C)+ 4H\l + u){B-C)- (4^^ + 24a'^) {B - D) - 2a\B" - D") 
3 

+^i/V2 + n){B' -C') + '^H\l + - C") - ^a''(i?' - D') 

oa a a 

+Ih^{2{1 + + 7Gi) + + + 5Gi') , (B.69) 

6 ba 

with D = du'C. Solving it for B- D, we get 

Among the nine bulk and six brane-boundary equations we started with, we have solved 
two bulk and one brane-boundary ones to determine Gi and B — D. Then we eliminate 
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Gi, G'l and B's from / / duE2, and E^ by (^l50|) and (pl70|) , and solve them for ^, C 
and G2 to get 



^ _ ^ 2Gi(p,) + 6(l + n)Gi(p,) + (18 + 3g,)Gi(p,) 
^ 6(2 + g,)g, 

n ^ - 3gj(l + ^)G'i(pj) - 6gjGi(pj) 

^ 6(2 + g,)g, ' ^ ' ^ 

• 2ij2(_6 + 2(-l + qj)u + (11 + gj)«2 + I2u^ + Sti^) 



6(2 + gj>2(2 + n 
^2(1 + -u)(i6 + qju{-A + 2u + Au"^ + n^)) 



2qj{2 + qj)u'^{2 + u) 

3(2 + 



2 r 

^2 = -E^o-7T2{(5i(2 + 9.)+3(4 + (?,)(l + n)2)Gi(^,,) 



+3u{l + u){2 + u)Gi{pj)Y (B.74) 

The last job is to verify the redundancy of the remaining four bulk equations and check if 
our solution satisfies the remaining five brane-boundary equations. It is easy to check 

duEj-Es = 0, (B.75) 



upon getting rid of Gi and Gi. As for E4, we can replace A by Z using ( B.36| ). The 



resulting equation has B, C, Gi and Z. Then use E2, E^ and E^ to write Z in terms 
oi B, C and Gi. The last step is to use ( [B.70 ) and ( [B.50| ), and we will see £'4 = 0. 



Similarly, we can show that Ei and Eq are redundant. Showing that our solution satisfies 
the brane-boundary equations is straightforward, using (pl^ ) to get the final answer. 
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